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ABSTRACT 

We study protoplanetary disc evolution assuming that angular momentum transport 
is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) 
turbulence in the hot inner regions. At radii of the order of 1 AU such discs develop 
a magnetically layered structure, with accretion occurring in an ionized surface layer 
overlying quiescent gas that is too cool to sustain MHD turbulence. We show that lay- 
ered discs are subject to a limit cycle instability, in which accretion onto the protostar 
occurs in ~ 10 4 yr bursts with M ~ lO~ 5 Af yr~ 1 , separated by quiescent intervals 
lasting ~ 10 5 yr where M « lO _8 M0yr _1 . Such bursts could lead to repeated episodes 
of strong mass outflow in Young Stellar Objects. The transition to this episodic mode 
of accretion occurs at an early epoch ((« 1 Myr), and the model therefore predicts 
that many young pre-main-sequence stars should have low rates of accretion through 
the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more 
massive than the minimum mass solar nebula, with most of the mass locked up in 
the quiescent layer of the disc at r ~ 1 AU. The predicted rate of low mass planetary 
migration is reduced at the outer edge of the layered disc, which could lead to an 
enhanced probability of giant planet formation at radii of 1 - 3 AU. 

Key words: accretion, accretion discs — MHD — stars: pre-main-sequence — stars: 
formation — solar system: formation — planets and satellites 



1 INTRODUCTION 

The structure and evolution of protoplanetary discs depend 
upon the rate at which gas can shed its angular momentum 
and thereby flow inwards. Two widely applicable physical 
mechanisms are known to lead to the required outward an- 
gular momentum transport. If the gas is coupled to a mag- 
netic field, instabilities that inevitably arise in differentially 
rotating discs (Balbus & Hawley 1991; Chandrasekhar 1961; 
Velikhov 1959) lead to turbulence and angular momentum 
transport (Stone et al. 1996; Brandenburg et al. 1995; for a 
review see e.g. Hawley & Balbus 1999). If the disc is massive 
enough, gravitational instability leads to additional trans- 
port (Toomre 1964; Laughlin & Bodenheimer 1994; Nelson 
et al. 1998; Pickett et al. 2000). 

Applying these findings to the construction of proto- 
planetary disc models leads to the structure shown schemat- 
ically in Fig. 1 (after Gammie 1996). In the inner disc, 
MHD turbulence transports angular momentum. However, 
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at larger radii of r ~ 1 AU, where the temperature is typi- 
cally a few hundred K, magnetic field instabilities are sup- 
pressed by the low ionization fraction (Matsumoto & Tajima 
1995; Gammie 1996; Gammie & Menou 1998; Livio 1999; 
Wardle 1999; Sano & Miyama 1999; Sano et al. 2000). This 
leads (Gammie 1996) to the formation of a layered disc 
structure, in which the gas near the disc midplane is cold, 
shielded from ionizing high energy radiation, and quiescent 
(non-turbulent). Turbulence and accretion occurs only in a 
thin surface layer that is ionized by cosmic rays. Moving still 
further outwards the entire thickness of the disc again be- 
come viscous, either at the radius where the surface density 
is small enough for cosmic rays to penetrate to the midplane, 
or where the onset of disc self-gravity provides an alternative 
non-magnetic source of angular momentum transport. 

The predictions of a static layered disc model for the ac- 
cretion rate and spectral energy distribution of T Tauri stars 
were discussed by Gammie (1996), and are broadly consis- 
tent with observations (e.g. with the accretion rate for Clas- 
sical T Tauri stars measured by Gullbring et al. 1998). In this 
paper we consider the evolution of the layered disc, which 
cannot be in a steady state (Gammie 1996, 1999; Stepinski 
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Figure 1. Illustration of the radial structure of the disc in the lay- 
ered disc model. At small radii, the central temperature exceeds 
Tcrit i the temperature above which MHD turbulence provides an 
efficient source of angular momentum transport. In the intermedi- 
ate, layered region, the central temperature is too low for the disc 
to support MHD turbulence. Accretion occurs through a surface 
layer which is kept ionized by cosmic rays, while the disc near the 
midplane is quiescent. At larger radii, the entire thickness of the 
disc is again viscous, with angular momentum transport being 
driven by self-gravity. 
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Table 1. Opacity regimes in ascending order of temperature, fit- 
ted by analytic functions of the form n = nop a T b .We have used 
fits provided by Bell & Lin (1994), as modified for low tempera- 
tures by Bell et al. (1997). The maximum temperature T max for 
each regime is quoted for a typical disc density of 10~ 9 gem -3 
(where specification of the density is necessary). 



1999), and examine the implications for the outflow history 
of young stars and for the predicted disc mass. The most 
significant changes to the disc structure occur at the radii of 
greatest interest for planet formation (Reyes-Ruiz & Stepin- 
ski 1995), and we discuss the implications for the migration 
of low mass planets, and for the eccentricity of massive plan- 
ets interacting with the disc. 



2 LAYERED PROTOPLANETARY DISC 
EVOLUTION 

2.1 Equations 

Describing the evolution of the surface density E(r, t) and 
midplane temperature T c (r, i) of a layered disc requires only 
minor modifications to the usual time-dependent equations 
for thin accretion discs. We denote the surface density of the 
'active' (viscous) disc by E a . If, 



(1) 



or, 

S <C 2 Slayer (2) 

then the disc is viscous throughout its thickness and E a = E. 
Otherwise only the surface layers are viscous and E a = 
2Ei ayor . The values of these parameters are determined by 
the requirement that the disc be sufficiently ionized to sup- 
port MHD turbulence (Gammie 1996). We adopt Tcrit = 
800 K, and Ei ayer = 10 2 gcm~ 2 . For a Keplerian disc, the 
angular velocity is = yj GM* / r 3 , where M* is the stellar 
mass. The surface density evolution is then described by, 



9E 
dt 



3d_ 

r dr 



1/2 



d 



,1/2^1 



+ E(r,t), 



(3) 



where v is the kinematic viscosity and E(r, t) is the rate of 
change of the surface density due to infall onto the disc. 

For the energy equation, we adopt a simplified form of 
that used by Cannizzo (1993), 



dTc 
dt 



2(Q+-Q_) 



dTc 
dr 



(4) 



Here c p is the disc specific heat, which for temperatures 
T c ~ 10 3 K is given by c p ~ 2.772.//U, where 1Z is the gas 
constant and /i = 2.3 is the mean molecular weight. Q + 
represents local heating due to viscous dissipation, given by, 



Q+ = (J) ^Efi 2 

if the entire disc is viscous and 
'9> 



Q + - (|) ,1 



'layer 



(5) 



(6) 



otherwise. For Q_, the local cooling rate, we assume that 
each annulus of the disc radiates as a blackbody at temper- 
ature T e , so that 



Q- = oTt, 



(7) 



where a is the Stefan-Boltzmann constant. Finally, we in 
elude an advective term in the energy equation, which de 
pends on the vertically averaged radial velocity, 

3 d 



V r 



i/Er 



l/2\ 



(8) 



ErV2 Qr 

and the radial temperature gradient. 
2.2 Vertical structure 

Completing the model requires specification of both the vis- 
cosity v and the vertical structure, which sets the relation 
between the central temperature T c and the surface tem- 
perature T e . We adopt the simplest, vertically averaged ap- 
proach, for which, 

3 



T 

± C 



(9) 



where r = (E/2)k is the optical depth for a given opacity 
n(p c , T c ). When an annulus makes the transition to the lay- 
ered state, we crudely account for this by replacing (E/2) 
in the expression for r by Ei aycr . Note that this means that 
we do not attempt to treat the vertical structure during the 
transition consistently. 

Analytic expressions for low temperature Rosseland 
mean opacities are given by Bell et al. (1997). The behaviour 
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of the disc depends primarily on the opacity near the tran- 
sition temperature T cr it, for which the fit is, 



k = 2x 10~ 2 T C 08 cm 2 g" 1 . 



(10) 



The full list of opacities used is quoted in Table 1. These fits 
have been taken from Bell & Lin (1994), with the modifica- 
tions for low temperatures quoted in Bell et al. (1997). 



2.3 Viscosity 

We adopt an alpha prescription (Shakura & Sunyaev 1973) 
for the viscosity, 



ft 



(11) 



where c s — y^TZTc/fi is the midplane sound speed and a a 
dimensionless parameter measuring the efficiency of angular 
momentum transport. Numerical simulations of MHD tur- 
bulence in discs suggest that a ~ 1CP 2 (Stone et al. 1996; 
Brandenburg et al. 1995). 



2.4 Self-gravity 

An immediate consequence of the layered disc structure is 
that the accretion rate through the surface layer is a fixed, 
increasing function of radius (Gammie 1996). Inevitably, 
therefore, material builds up in the quiescent layer over 
time, until the surface density becomes high enough that 
self-gravity of the disc becomes important. The condition 
for this is that the Toomre Q parameter (Toomre 1964), 
given by, 



Q = 



c s ft 



(12) 



is less than some critical value Qcrit, which is of the order 
of unity. Here, we adopt Q cr it = 2. When Q < Q CI it, non- 
axisymmetric instabilities set in and act to transport angu- 
lar momentum outwards. In doing so, the surface density 
evolves in such a way that the stability of the disc generally 
increases, thereby self-regulating the violence of the insta- 
bility. 

In a local approximation, the efficiency of angular mo- 
mentum transport due to self-gravity can be determined 
by the requirement that viscous heating balances radiative 
losses (Gammie 1999). This approach, in which the cool- 
ing time directly determines a, is most appropriate for discs 
whose thermal balance is determined by viscous heating. 
Even in this limit, global effects may be important (Balbus 
& Papaloizou 1999). A global redistribution of surface den- 
sity, driven by self-gravity, can also occur even if the disc 
temperature is set externally (for example by irradiation). 
This is the limit that has been studied extensively using 
isothermal or polytropic simulations (e.g. Laughlin & Bo- 
denheimer 1994). For protoplanetary discs, models suggest 
that both viscous heating and irradiation can be important 
at different epochs and radii (Bell et al. 1997). 

Given the above complications, it is evident that at- 
tempts to include the effects of self-gravity into one dimen- 
sional disc models can only be approximate. We adopt the 
form suggested by Lin & Pringle (1987, 1990), and use a 
variant of equation (Till) along with an effective gravitational 



'a', i/ grav = a g ravC 2 /f2, with 

Q grav = 0.01 - 1 j if Q < Qcrit 

Ctgrav = otherwise. (13) 

Operationally, we solve equations (^) and (Q) using an ex- 
plicit finite difference method, and include the self-gravity 
as an additional sweep over the grid using f grav and the full 
surface density £ in the diffusive part of equation (^|), and 
as an additional heating term in the energy equation. 



3 RESULTS 

For our numerical calculations we adopt M* = Mq , an inner 
disc radius r ln = 3.5 x 10 11 cm (5 Rq), and an outer disc 
radius r out = 6.0 x 10 14 cm (40 AU). The computational 
grid of 120 radial mesh points is uniform in a scaled radial 
variable X <x r 1 ' 2 . A zero torque boundary condition is im- 
posed at r- m , while at r out we prevent outflow by setting the 
radial velocity, v r , to zero. When mass is added to the disc 
(representing infall of further material from the molecular 
cloud), we take E to be a Gaussian centered at 10 AU with 
a width of 1 AU. This is exterior to the layered section of 
the disc that we are interested in studying, so the details 
of how mass is added are unimportant here. We add mass 
assuming that it has the same specific angular momentum 
and temperature as the local disc material. 



3.1 Steady infall models 

To demonstrate the evolution of the layered disc model, we 
first take the infall rate onto the outer disc Mi n f a n to be 
a constant. We take a — 10~ 2 , as suggested by numerical 
simulations of MHD disc turbulence, and allow the disc to 
evolve from an initially low mass (10 -2 Mq) until either a 
steady state or a limit cycle has developed. 

Figure 2 shows the accretion rate onto the star for a 
range of mass infall rates between 5 x 10~ 7 Moyr -1 and 
3 x 10~ 6 M0yr _1 . For the highest infall rate, the disc is able 
to transport gas steadily onto the star. The inner region, in 
which T c > T cr it, is able to be supplied by an outer disc in 
which Q < Qcrit and the angular momentum is transported 
by self- gravity. Q is typically around 1.5 in this region, with 
the transition occurring at a radius of ~ 3 AU. 

For lower infall rates such steady accretion is not pos- 
sible. After an initial phase in which the disc mass increases 
steadily, a limit cycle is obtained in which outbursts, with 
a peak accretion rate M « 10 -5 Mqvt -1 , alternate with 
quiescent intervals where M « 10 -8 Mqyt' 1 . Both the du- 
ration of the outbursts, and the recurrence time, vary with 
the rate at which gas is supplied to the outer disc. However, 
to order of magnitude, for Mi n f a n < 10 -6 Moyr -1 we find 
^outburst ~ 10 4 yr, and t ICCUT ~ 10 5 yr. 

The mechanism for these outbursts is shown in Figure 
|§| which plots the central temperature during the quiescent 
interval leading up to an outburst. Following the end of one 
outburst, there is an extended layered region of the disc be- 
tween ~ 0.3 AU and ~ 3 AU. This region is stable against 
gravitational instability, and supports only the small qui- 
escent accretion rate. Gas flowing inwards from the more 
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Figure 2. Protostellar accretion for models with a constant 
rate of infall onto the outer disc. From top down, the pan- 
els show models with M; n f a n = 3 X 10~ 6 MQyr -1 , Mj n fall = 
1.5X1CT 6 M Q yr _1 , and Minfall = 5xl0~ 7 M Q yr _1 respectively. 
All models are shown after initial transients have decayed. Steady 
accretion occurs for sufficiently high accretion rates through the 
outer disc. At lower accretion rates, the mass flow onto the star 
is strongly time-dependent. 



Figure 3. Evolution of the disc central temperature during qui- 
escence. An outburst is triggered when the central temperature 
in the layered region first exceeds T cr ; t (dashed line) , the temper- 
ature at which the disc is well enough ionized to support MHD 
turbulence. The curves show temperature profiles plotted at in- 
tervals of 10 4 yr, with the earliest timeslice having the lowest 
temperature at 2 AU. The final timeslice, shown as the upper, 
bold curve, immediately precedes the start of an outburst. 



active self-gravitating region of the disc thus accumulates in 
the quiescent region, and the transition radius where the disc 
becomes self-gravitating moves inward. At the same time, 
the greater dissipation rate in the self-gravitating disc cre- 
ates a local peak in the disc temperature, which also moves 
inwards. A new outburst is triggered when this peak first 
exceeds T cr it, thereby satisfying the condition for MHD tur- 
bulence to be restarted. Once the outburst begins thermal 
energy is rapidly advected inwards, triggering a transition 
of the entire inner disc into a hot turbulent state with high 
accretion rate. This continues until the reservoir of accu- 
mulated mass has been depleted, whereupon a cooling wave 
sweeps through the disc and returns it to the quiescent, lay- 
ered, state. 

The derived outburst durations are consistent with the 
expected viscous timescale in protoplanetary discs at radii 
of a few AU. For an outburst triggered at a radius of ap- 
proximately 2 AU, the viscous timescale once the disc is in 
outburst, given by 

^i7=^y * (i4) 

where h is the disc scale height, is approximately, 

<--^rww v <i5) 

This is only a rough estimate, but it is of the correct order 
of magnitude to match the numerical model. 

The involvement of disc self-gravity means that rela- 



tively large disc masses are required for the limit cycle to 
operate. In our models, Mdi sc fluctuates between 0.2 and 
0.3 Mq. During an outburst, only a small fraction of the 
disc mass — around 20% — is accreted. In more realistic 
evolutionary models, in which the mass infall rate declines 
with time, the disc will therefore be able to produce a hand- 
ful of outbursts, with a few hundredths of a solar mass of 
gas accreted during each, before settling into a permanent 
quiescent state. 

3.2 Relation to FU Orionis events 

The most striking variability observed in pre-main-sequence 
stars occurs in FU Orionis events (e.g. Kenyon 1995; Hart- 
mann & Kenyon 1996) , which are large amplitude outbursts 
in the system luminosity originating in the accretion disc. 
The statistics of these events are subject to considerable un- 
certainty, but the peak accretion rate during outbursts is of 
the order of 10 -4 M Q yr _1 , while the duration is generally 
thought to be around 10 2 yr. 

The most popular explanation for FU Orionis events is 
in terms of a disc thermal instability (e.g. Bell & Lin 1994; 
Kley & Lin 1999, and references therein), akin to that used 
to model dwarf novae (e.g. Cannizzo 1993). In the protostel- 
lar case, a thermal instability can only operate at extremely 
small disc radii, typically at less than 0.1 AU. Matching the 
observed timescales then requires small values of a (10~ 3 
to 10~ 4 ), with correspondingly large values for the disc sur- 
face density. As noted by Gammie (1999), it would be at- 
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tractive to dispense with the need for a thermal instabil- 
ity by invoking an alternative limit cycle of a layered disc. 
This might be possible if the transition to the outburst state 
was triggered at the extreme inner edge of the layered re- 
gion. Whether this happens depends on where in the quies- 
cent layer mass preferentially accumulates, and the outcome 
is therefore sensitive to the detailed, and rather uncertain, 
physics of the layered disc. In our model, in which matter is 
being added to the layered region primarily at large radius, 
the triggering occurs at the outer edge of the layered region, 
and the timescales are inconsistent with those of FU On- 
onis events. We note, however, that the accretion rates of 
(1 - 10) x Hn 5 M yr _1 that we obtain during our outbursts 
fall into the thermally unstable regime. The long outbursts 
obtained here could then feed shorter FU Orionis event of 
the sort calculated by Bell & Lin (1994). 

3.3 Protostellar accretion history 

Continuous replenishment of the disc mass from infall is not 
a realistic model for protostellar accretion. To explore how 
the accretion rate onto the star may evolve with time, we 
adopt a simple model for the infall of gas onto the disc. The 
details are model-dependent (see e.g. Larson 1969; Shu 1977; 
Basu 1998), but at early times the infall rate is expected to 
be of the order of Minfaii ~ <? S /G, where c s is the sound 
speed in the collapsing cloud. For cloud temperatures T ~ 
10 K, this implies an infall rate of ~ 10 -5 M^yr -1 . We take 
an initial infall rate of 2 x 10~ 5 M©yr — 1 , and assume that 
this declines exponentially on the free fall timescale, tg = 
(37r/32Gp c i ud) 1/ ' 2 , where p c i ou d is the cloud density. We take 
ts = 10 yr. Other input parameters are as described earlier. 

With these parameters, Fig. ^ shows the accretion rate 
for the layered disc model. For comparison, we also show the 
results from a viscous disc model which has identical param- 
eters, but which has no layered region. This control model 
would be appropriate, for example, if angular momentum 
transport was driven by a purely hydrodynamic mechanism 
whose efficiency was independent of the disc temperature. 
We compute the control model by using the same code and 
infall history, but with T cr it set to zero. 

Initially, when the infall rate is higher than around 
2 x 10 -6 MqJT -1 , both models support an accretion rate 
onto the star which tracks that with which matter is be- 
ing added to the disc. Subsequently, accretion through the 
layered disc occurs in outbursts with properties similar to 
those described earlier, but with increasing intervals between 
events. For this model, the last outburst before the disc be- 
came permanently quiescent occurred after 2 Myr. The ac- 
cretion rate onto the star in the control model, on the other 
hand, declines smoothly with time. Similar results are ob- 
tained using other functional forms for the mass infall rate 
as a function of time. 

The timescales of variability predicted by the model are 
much longer than any direct observational record. Indirect 
evidence for long timescale variability of young stars, how- 
ever, is provided by studies of protostellar outflows, since the 
rate of mass outflow via jets is widely believed to track the 
accretion rate. For example, observations by Reipurth, Bally 
& Devine (1997) suggest that large amplitude variations in 
the outflow rate occur not only on the 10 2 yr timescales 
characteristic of FU Orionis events, but also on much longer 




Figure 4. Accretion rate and disc mass as a function of time. 
The solid curves are for layered disc models in which the mass 
infall rate onto the outer disc declines exponentially. For compar- 
ison with these, the dashed curves show the evolution of a model 
without a layered structure. This model has identical parame- 
ters, except that the disc is assumed to be viscous at all radii and 
central temperatures. 



timescales of 10 4 yr or greater. These longer timescales are 
characteristic of processes occurring at larger disc radii than 
thermal instabilities, and are consistent with instabilities, of 
the kind discussed here, originating in the layered disc re- 
gion. 

The outburst accretion rate in the layered disc model is 
several orders of magnitude higher than that in the equiva- 
lent fully viscous disc model at t ~ 10 6 yr. However, the time 
averaged accretion rate is substantially lower, so that at late 
times the disc mass in the layered model substantially ex- 
ceeds that in the viscous model. We find that Afdisc ~ 0.1 Mq 
after 1-2 Myr. By this time the mass of the fully viscous 
disc has dropped to only 10 -2 Mq. As shown in Fig. 5, most 
of this extra mass is tied up in the quiescent layer at radii 
of the order of 1 AU, where the surface density exceeds the 
minimum mass solar nebula estimate (Hayashi, Nakazawa 
& Nakagawa 1985) by about two orders of magnitude. 

3.4 Compatibility with observations of 
protoplanetary discs 

A disc mass of a few tenths of a solar mass at early times is 
consistent with the upper end of the disc mass distribution 
inferred from mm-wavelength observations (Beckwith et al. 
1990; Osterloh & Beckwith 1995). However, the layered disc 
model predicts that a substantial mass - at least a few hun- 
dredths of a solar mass - remains in the disc throughout 
most of the typical Classical T Tauri disc lifetime (Strom 
1995). Although recent measurements of the mass of H2 in 
debris discs indicate that several Jupiter masses of gas can 
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r / AU 

Figure 5. Cumulative disc mass as a function of (upper panel) 
surface density S and (lower panel) radius. The solid curve shows 
the layered disc model at t = 1.5 X 10 yr. For comparison, the 
dashed curve shows results for the fully viscous model at the 
same time. The vertical dashed line in the upper panel shows 
the approximate surface density below which the disc would be 
optically thin at A = 1 mm. Essentially all the additional mass 
present in the layered disc is locked up in highly optically thick 
regions at small radii. The amount of mass at r > 10 AU is 
comparable in both the layered and the fully viscous model. 



persist beyond 10 Myr (Thi et al. 2001), our predicted disc 
masses are in apparent contradiction with mm observations 
of many systems with disc masses estimated at 10~ 3 Mq or 
smaller. 

To check whether this is a serious problem for the lay- 
ered disc model, we plot in Fig. 5 the cumulative distribu- 
tion of disc mass as a function of surface density. Most of 
the 'additional' mass in the layered disc, compared to the 
equivalent fully viscous control model, is locked up in the 
quiescent layer at small radii and high surface density. This 
material is highly optically thick, and thus would not con- 
tribute to the observed emission at mm wavelengths. For 
a range of dust models, the opacity of dust at A = 1 mm 
is thought (Menshchikov & Henning 1997; Beckwith et al. 
1990) to be in the range 0.3 cm 2 g _1 < k < 10 cm 2 g~ x (note 
that this is per gram of dust). For a gas to dust ratio of 
100, this implies that an optical depth of unity at 1 mm is 
attained for S ~ 10 2 gem -2 . At this surface density thresh- 
old, the masses of optically thin gas in the layered and non- 
layered disc models are very similar. The masses of proto- 
planetary discs inferred from mm- wavelength measurements 
could thus substantially underestimate the true disc masses. 

The accretion rates inferred for Classical T Tauri stars 
can also provide a test. In the simplest version of the layered 
disc model, the accretion rate during quiescence does not 
change over the lifetime of the protostellar phase. This rate is 
consistent with the order of magnitude inferred for accreting 



Classical T Tauri stars (Gullbring et al. 1998). There have 
also been suggestions that the observationally inferred accre- 
tion rate declines with increasing stellar age (Hartmann et 
al. 1998). Although currently the scatter in measured accre- 
tion rates at a given assumed age is very large, and the age 
estimates themselves are subject to substantial uncertainties 
(Tout, Livio & Bonnell 1999), this observation poses diffi- 
culties for the simplest layered disc model (Stepinski 1998). 
We note, however, that changes in the quiescent accretion 
rate with time would be expected in more complete models, 
for example if the opacity from dust decreased with time 
due to grain growth. Including heating of the disc from the 
changing stellar radiation would also lead to a decline in the 
quiescent accretion rate. 



4 PLANET FORMATION 

The environment which a layered disc presents for planet 
formation differs in several respects from both previous vis- 
cous disc models (e.g. Lin & Pringle 1990), and from static 
models such as the minimum mass solar nebula. At radii be- 
tween a few tenths of an AU, and several AU, the gas near 
the disc midplane is almost always quiescent, unless the set- 
tling of cm sized particles itself drives additional instabilities 
(Goldreich & Ward 1973; Cuzzi, Dobrovolskis & Champney 
1993). The low viscosity of the disc in this region would 
reduce the rate at which massive planets migrate inwards, 
and may allow planet-disc interactions to excite significant 
eccentricity (Papaloizou, Nelson & Masset 2001). 

Other differences include the disc at large radii remain- 
ing mildly unstable to gravitational instability for a rela- 
tively long time - around a Myr - because the reduced time- 
averaged accretion rate leads to a larger disc mass at late 
times when compared to viscous disc models. Outbursts, 
which persist for a substantial fraction of the disc lifetime, 
mean that much of the inner disc is expected to be subject 
to rapid heating and cooling episodes. However, we find no 
evidence that these cooling waves drive the disc into a state 
where gravitational collapse (Q < 1) would occur. 

The survival of solid bodies in protoplanetary discs is 
limited by the rate at which they migrate inwards relative 
to the gas. Migration can be rapid both for cm-sized bodies 
(e.g. Godon & Livio 1999, and references therein), and is 
particularly problematic for low mass planets, for which mi- 
gration occurs due to the influence of gravitational torques 
(Goldreich & Tremaine 1979). In standard disc models, the 
migration timescale for Earth mass planets at a few AU 
can be as short as 10 4 — 10 5 yr, which leaves little time to 
assemble the cores of giant planets before the putative build- 
ing blocks are consumed by the star. The rate of migration 
depends only very weakly upon the surface density profile 
(Ward 1997), but is sensitive to the gradient of the central 
temperature, and can be halted if there exists a region of the 
disc where T c oc r. During the quiescent phase, the models 
we have discussed here possess such a non-monotonic cen- 
tral temperature profile at radii of 1-3 AU (Fig. 3), leading 
to the possibility of an enhanced probability of giant planet 
formation at those radii (Papaloizou & Terquem 1999). Cur- 
rent radial velocity surveys are sensitive to massive planets 
at radii r < 3 AU, i.e. within this region. If migration is 
as rapid as currently suspected, a consequence of the lay- 
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ered disc model would probably be fewer planets at larger 
radii than would be expected from an extrapolation from 
the current data. 



5 SUMMARY 

In this paper, we have presented models for the evolution of 
magnetically layered protoplanetary discs. Given our present 
understanding of angular momentum transport in discs, this 
model represents the best guess for the structure of proto- 
planetary discs at radii of the order of 1 AU, where the gas is 
cool and poorly coupled to the magnetic field. The evolution- 
ary models presented here suggest a number of important 
differences with non-layered viscous disc models. 

(i) The disc cannot be in a steady state for outer disc 
accretion rates in the range 1CP 8 Moyr -1 < M < 2 x 
1CP 6 Moyr -1 . A limit cycle is obtained, in which heating 
when the layered region of the disc becomes self-gravitating 
periodically restarts MHD turbulence, leading to outbursts 
of accretion onto the star. All calculations to date of layered 
discs obtain strongly episodic accretion, despite variations 
in the physical processes included in the models (Gammie 
1999; Stepinski 1999). This is the most robust prediction of 
the model. 

(ii) The duration of outbursts depends on where in the 
quiescent disc they are triggered. In our model, the outbursts 
are long, with duration ~ 10 4 yr. They could drive repeated 
strong episodes of mass outflow from the inner disc, resulting 
in 'pulsing' of the observed jets. 

(iii) The time-averaged accretion rate is reduced due to 
the bottleneck created by the layered region of the disc. As a 
result, the disc mass at late times is large, typically « O.1M 
at 1-2 Myr. Most of this mass is locked up in the quiescent 
layer of the disc at small radii. 

(iv) The central temperature is strongly modified, and 
often increases with radius, near the transition between the 
layered disc and the outer self-gravitating region. This could 
slow the otherwise rapid rate of low mass planetary migra- 
tion at radii r ~ 1 — 2 AU. 

(v) A low viscosity in the layered region of the disc af- 
fects the migration rate and eccentricity evolution of mas- 
sive planets at the radii currently probed by radial velocity 
surveys. 

Current observations provide only very limited constraints 
on the properties of protoplanetary discs at the radii of 
greatest interest for planet formation. As a result, there re- 
mains substantial uncertainty in our knowledge of how these 
discs evolve. Models, such as this one, that attempt to in- 
clude more realistic disc physics, can lead to very different 
environments for planet formation and migration than the 
highly simplified models usually considered. The most ob- 
vious observational predictions of the model are that some 
very young stars should be accreting at the low 'quiescent' 
rate of the order of 10 -8 M0yr _1 , and that high accretion 
rate outbursts should continue, albeit with lesser frequency, 
during a substantial fraction of the Classical T Tauri phase. 
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